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Robert  W.  Courter 
ABSTRACT 

A  simple  time-dependent  model  of  the  internal 
ballistics  of  a  two-stage  light  gas  gun  is  briefly 
described.  The  gasdynamic  flow  field  between  the  driving 
piston  and  the  driven  sabot  is  singled  out  for  special 
study.  This  flow  regime  is  idealised  as  a  quasi-one- 
dimensional  flow  of  an  inviscid#  compressible#  real  gas 
with  arbitrarily  moving  longitudinal  boundaries.  A  method 
of  numerical  solution  of  the  governing  partial  differential 
equations  in  suggested.  in  addition#  a  resoning  technique 
for  describing  the  time-dependent  computational  grid  is 
presented. 
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I 


INTRODUCTION 


Ballistic  testing  in  an  instrumented  £ree-£light 

* 

range  is  an  important  part  of  the  broad  £ield  of 
experimental  aerodynamics.  Launcher  technology,  in  turn, 
is  a  vital  part  of  ballistic  testing.  The  flight  speed, 
initial  acceleration  history  and  exit  flow  environment  are 
all  dependent  on  launcher  design  parameters.  it  is 
important  in  ballistic  testing,  therefore,  to  have  the 
capability  of  predicting  launcher  performance. 

The  two-stage  light  gas  gun  is  one  example  of  a  large 
class  of  compressed  gas  launchers.  Such  a  gun  is  to  be 
used  at  the  Air  Force  Armament  Laboratory  as  a  primary 
launcher  for  a  variety  of  free-flight  model  configurations. 
A  satisfactory  performance  prediction  program  for  this 
launcher  is  required.  One  of  the  objectives  of  the 
research  pursued  during  this  summer  research  program  is  to 
develop  such  a  procedure.  A  second  objective  is  to 
investigate  the  complex  gasdynamic  processes  which  occur 
within  the  launcher  during  a  firing  cycle. 

A  literature  survey  of  DOD,  NTIS  and  NASA  data  bases 
revealed  only  sparse  information  on  important  aspects  of 
launcher  performance  such  as  piston  extrusion  behavior, 
diaphragm  rupture  dynamics,  launcher  experimental  data  and 
internal  flow  field  descriptions.  Gun  analysis  programs 
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for  other  facilities  have  been  developed  (see  References  1 
through  4)  which  yield  adequate  predictions  of  internal 
pressure  and  muzzle  velocity.  However,  these  codes  are 
either  outdated  or  unavailable,  or  both.  Also,  they 
apparently  do  not  provide  performance  without  considerable 
computation  expense,  nor  do  they  provide  detailed 
information  on  the  internal  flow  field. 

This  fundamental  data  deficiency  suggests  that 
the  formulation  of  a  very  simple  performance  prediction 
model  would  be  beneficial  and  that  detailed  studies  of 
the  internal  flows  would  add  to  the  understanding  of 
launcher  performance. 

In  the  present  study,  a  simple,  time-dependent 
volumetric  model  of  the  gas  gun  operation  is  incorporated 
into  a  computer  program  in  a  modular  form  which  will  permit 
investigation  of  various  aspects  of  the  launcher  cycle  in 
an  independent  manner.  In  addition,  more  detailed  models 
of  the  basic  flow  kinematics  and  thermodynamics  are  derived 
which  have  promise  in  revealing  the  gas  property 
distributions  during  gun  operation. 

II.  OBJECTIVES 

The  principal  objective  of  the  present  work  is  to 

develop  a  performance  prediction  model  for  a  class  of 

ballistic  launchers  which  includes  the  two-stage  light  gas 

gun.  A  typical  configuration  of  such  a  launcher  is  shown 

in  Figure  1.  It  is  intended  not  only  to  develop  a  model 
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which  can  be  useful  immediately  in  assessing  launcher 
performance  as  a  function  of  launcher  and  projectile 
geometric  characteristics  and  propellant  charge  weight,  but 
also  to  investigate  means  of  providing  a  more  detailed 
description  of  the  gasdynamics  within  the  gun. 

The  author  has  assisted  his  graduate  student,  Mr. 
Raymond  Patin,  in  the  formulation  of  a  baseline  launcher 
performance  model  which  can  be  used  in  the  initial 
predictions.  Parametric  studies  and  comparisons  with 
experimental  data  indicate  that  this  model  is  reasonable. 
Details  of  the  basic  model  are  presented  in  Reference  5. 
However,  essential  model  features  are  outlined  in  this 
report  to  provide  background  for  the  discussion  of  flow 
analysis  methods  to  be  discussed  here.  In  the  present 
report  details  are  given  of  the  analytical  procedures 
necessary  to  determine  the  quasi-one-dimensional 
distributions  of  thermodynamic  and  kinematic  properties  of 
a  flow  field  confined  by  two  boundaries  moving  at  different 
speeds. 

III.  THE  BASELINE  LAUNCHER  PBRFORMAtlCE  MODEL 

In  Figure  1  are  shown  the  physical  details  of  a 
typical  two-stage  light  gas  gun.  The  gun  operation 
sequence  begins  with  the  ignition  of  a  single-perforation 
propellant  charge  in  the  combustion  chamber.  The  pressure 
rise  created  by  gas  generation  during  combustion  forces 

the  piston  to  move  toward  the  diaphragm,  compressing  the 
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helium  gas  in  the  process*  At  some  point  during  this 
compression  process,  the  scored  metal  diaphragm  ruptures, 
exposing  the  base  of  the  sabot  to  high-pressure  gas.  The 
sabot  is  then  propelled  at  high  speed  through  the  launch 
tube* 

In  Reference  S  a  simple  model  of  this  launching 
process  is  described  in  detail.  Only  an  outline  of  this 
performance  model  will  be  given  here.  The  time- 
dependent  equations  which  are  integrated  to  provide  a 
simulation  of  the  launcher  operation  are  listed  in 
functional  form  in  Table  1.  The  pressure  variation  in  the 
combustion  chamber  is  described  by  a  differential 
equation  derived  by  0.  K.  Heiney  in  Reference  6.  The 
Kinematics  of  the  piston  and  sabot  are  determined 
from  the  conventional  translational  momentum  equations. 
The  gas  properties  within  the  pump  tube  and  launch 
tube  chambers  are  determined  by  isentropic  pressure/volume 
relationships. 

This  basic  prediction  program  appears  to  describe  the 
trends  of  launcher  operation  correctly.  However,  to  date, 
insufficient  comparisons  of  the  model  predictions  with 
experimental  data  have  been  made  to  declare  it  a  valid  and 
versatile  model.  One  of  the  weaknesses  of  this  baseline 
prediction  method  is  that  it  cannot  account  for  the  details 
of  the  gasdynamics  within  the  pump  tube  of  the  launcher. 
In  the  following  sections  of  this  report  methods  for 
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detailed  analysis  of  this  flow  regime  will  be  discussed. 

IV.  GASDYNAMICS  ANALYSIS  METHODS 

The  flow  field  in  the  pump  tube  section  of  a  two- 
stage  light  gas  gun  can  be  idealized  as  a  quasi-one 
dimensional  flow  in  a  tube  of  variable  area  with 
longitudinal  boundaries  which  are  moving  at  different 
rates.  Of  particular  interest  in  this  class  of  flows  are 
the  characteristics  in  any  regions  of  variable  area  and 
in  the  vicinity  of  the  moving  boundaries.  In  the  case  of 
the  light  gas  gun  this  latter  aspect  is  further  complicated 
by  the  fact  that  one  of  the  boundaries  (the  piston)  may 
intrude  into  a  region  of  decreasing  area  (extrude) , 
resulting  in  a  surface  velocity  which  is  not  the  same  as 
the  velocity  of  the  center  of  mass  of  the  piston.  In  any 
flow  involving  high  pressures  and  temperatures  the  effects 
of  departures  of  gas  behavior  from  the  ideal  gas  law  model 
must  also  be  evaluated. 

In  this  section  of  the  report  an  anlytical  model  for 
the  flow  of  a  real  gas  in  a  variable  area  tube  with  moving 
longitudinal  boundaries  is  presented.  No  attempt  is  made 
to  solve  the  equations  at  this  time,  but  a  method  of 
solution  is  presented  and  discussed. 

The  quasi-one-  dimensional  equations  of  flow  for  an 
inviscid,  compressible  fluid  can  be  written  in  conservative 
form  as  follows: 


CONSERVATION  OF  MASS 


o 


(1) 


.  m  4A  m 

T  A  dx 


CONSERVATION  OF  MOMENTUM 


CONSERVATION  OF  ENERGY 


(2) 

(3) 


where  f  is  the  mass  density  (mass  per  unit  volume),  m  is 
the  translational  momentum  per  unit  volume,  e  is  the 
total  energy  per  unit  volume,  p  is  the  thermodynamic 
pressure,  X  is  the  axial  coordinate,  t  is  the  time  and 
A  is  the  cross-sectional  area  of  the  tube.  The  total 
energy,  e  ,  is  given  by 


where  i  is  the  internal  energy  per  unit  mass  of 
the  gas.  The  gas  in  the  confined  chamber  is  assumed  to  be 
helium.  An  appropriate  real  gas  model  for  helium  at  high 
pressure  and  temperature,  given  by  Harrison  in  Reference  7, 
is  proposed  for  use  here.  This  model  is  a  fourth  order 
virial  formulation  given  as  follows: 


P  =  jRt[  i  +  B(r)f  +  cfrtf1  +  o60f3] 


(5) 


where  R  is  the  engineering  gas  constant  and  B (T) ,  C(T)  and 
D (T)  are  empirical  functions  of  the  temperature. 

Equations  (1)  through  (5)  govern  the  quasi-one- 
dimensional,  unsteady,  inviscid  flow  of  a  real  gas 
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(helium).  The  boundary  conditions  for  the  flow  of  interest 
here  are  as  follows: 

At  X  *  Xp  ,  the  boundary  moves  with  velocity  xp 
At  X  *  L  +  Xs  ,  the  boundary  moves  with  velocity  Xs 

where  L  is  the  initial  separation  of  piston  and 

sabot. 

The  solution  of  the  partial  differential 
equations  and  boundary  conditions  which  govern  this  flow 
situation  must  be  effected  by  numerical  methods.  In  the 
following  section  an  appropriate  method  of  numerical 
integration  for  the  equations  is  given.  Also,  a  procedure 
for  dealing  with  computational  grid  dynamics  by  resoning  is 
presented. 

IV.  NUMERICAL  METHODS 

In  1969  MacCormack  published  a  second-order  accurate 
numerical  technique  for  the  solution  of  a  class  of  partial 
differential  equations  (see  Reference  8) .  The  method,  used 
successfully  by  the  author  in  a  multidimensional  gasdynamic 
analysis  (Reference  9) ,  is  applicable  to  the  set  of 
hyperbolic  partial  differential  equations  which  govern  the 
present  flow  study.  In  general  terms,  the  two-level 
technique  can  be  written: 

» 

Governing  Equation 

+  3X  +  G  *  O  (6) 
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where  U  is  the  conserved  function  vector  and  F  and  G 


are  the  flux  vectors. 


Numerical  Solution 


Predictor 


F»  “  F;-, 


=  4  -  At  +  fi 


Corrector 


_ _  _  _<H>  _ _  _ 

U-*'=  z  [u"+  u;-  At  ( 'iT7J  *  O]  I’) 

where  n  is  the  time  index  and  j  is  the  space  index. 

This  explicit  method  is  very  simple  to  set  up  for  the 

particular  set  of  equations  to  be  solved  in  the  present 

case.  These  equations  are: 


CONSERVATION  OF  MASS 


Corrector 


Predictor  n+i  *  /rhj-Rfy.,  in*  d A  \ 

?j  *  fj  ~  \  ^  A  j  dtf  j  / 


CONSERVATION  OF  MOMENTUM 


CION  OF  ENERGY 

’» .  j-«na*faBaC.  itfS  j 


Corrector 


CONSERVATION  OF  ENERGY 


pt,dict6f 


C0,Trt0^(e>7LAt{ff^C-ffrf4  (1” 

The  particular  solution  of  interest  here  involves  a 


flow  field  between  boundaries  which  are  moving  at 


arbitrary,  different  speeds.  This  presents  some 
computational  difficulties.  The  computational  grid  which 
must  be  established  to  effect  the  numerical  solution  of 
Equations  (8)  through  (10)  must,  necessarily,  be  in  motion.' 
in  order  to  produce  a  result  which  is  accurate  the  grid 
must  deform  and  move  simultaneously.  This  effect  can 
be  achieved  by  the  formulation  of  an  "adaptive"  grid.  One 
method  of  constructing  such  a  grid  for  a  one-dimensional 
problem  is  by  the  rezoning  technique  which  will  be  used  in 
the  present  solution.  "Rezoning"  implies  that  at  the 

end  of  each  computation  step,  the  computational  grid  is 
adjusted  to  conform  to  the  new  boundary  conditions.  Such 
adjustments  may  involve  tedious  calculations  at  each  time 
step  to  prevent  the  transmission  of  numerically-induced 
flow  perturbations,  but  they  ensure  a  solution  with 
appropriate  resolution.  Consider  the  sketches  of 

Figure  2.  In  part  (a)  is  shown  the  grid  structure  which 
was  used  to  advance  the  solution  from  time  step  "n"  to  time 
step  "n+l".  In  that  time  interval  the  computational 
boundaries  were  altered  by  the  boundary  motion.  In  part 
(b)  is  shown  the  adjusted  structure  (the  rezoned  structure) 
for  time  step,  "n+l".  The  equations  which  describe  the 
grid  distribution  and  the  corresponding  values  of  the  state 

variables  can  be  given  as  follows  (refer  to  Figure  2) : 

# 

COMPUTATION  DOMAIN 

Lm,=  L  -  (U-Xs)  (ll» 
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CELL  SIZE 


mi  n 

•4  * 

CELL  BOUNDARY  LOCATION 


AX, 


X4 


mi 


n+i 


,  v  mi  m« 

/*£  -**  \ 

-  V  N  / 

(12) 

i-' 

z 

It*  i 

(13) 

STATE  VARIABLE  ALLOCATION  (Typical  Internal  Cell) 
,  Hi!  /*«  .mi  ,  mi.  mi  \/  t\ 

4>j  =  Uj  <t>i-  +  b;  <t>i  ;Axj 


where 


and 


n+i 


<>i  =  *r+  z  (Axr-axo 

if  -  £<  -M£  ax?1-  xPnw 

4f 


(14) 


Ha*  tt»i 

is  the  state  variable  value  at  the 
grid  point  for  time,  t<*'  . 


4> 


fVH  is  the  state  variable  value  at  the 
grid  point  for  time,  t1*1  . 


new 


old 


These  equations  are  a  typical  set.  A  different  set  is 
required  for  each  different  pattern  of  boundary  motion. 
These  equations  have  been  established  based  on  a  uniform 
grid.  Grid  clustering,  which  may  be  necessary  in  some 
cases,  will  be  discussed  in  the  next  section  of  the 
report. 


It  is  important  to  establish  the  sequence  of  solution 
of  the  above  set  of  equations.  The  computational  steps  to 
be  employed  at  each  grid  point  in  advancing  the  solution  • 
from  one  time  step  to  the  next  are: 

1.  Determine  the  computational  domain  and  establish  the 
grid  size. 

2.  Determine  the  computational  time  step  by  invoking 
the  Courant  condition. 

3.  Advance  the  state  variables  to  the  next  time  by 
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application  of  the  MacCormack  algorithm. 

4.  .Solve  the  dynamic  equations  of  the  boundaries  for 

the  given  time  step.  (For  the  case  of  the  light 
gas  gun  this  involves  solving  the  kinematic 
equations  which  are  listed  in  Table  1.) 

5.  Perform  the  rezoning  operations. 

6.  Continue  with  the  solution. 

V.  ADVANCED  CONSIDERATIONS 

The  foregoing  computational  methods  have  been 

functional,  but  not  elegant.  There  are  other  methods  which 

will  result  in  a  more  efficient  code  and,  consequently, 

give  better  and  quicker  results.  The  MacCormack 

integration  algorithm,  like  all  second  order  systems,  has  a 

serious  deficiency  in  regions  of  rapid  change  of  the  state 

variables  (such  as  near  shock  waves)  ...  it  becomes 

divergent.  The  use  of  artificial  diffusive  terms  in  the 

equations  prevents  the  divergence,  but  it  results  in  an 

oscillatory  solution  and  tends  to  smear  the  effective 

discontinuities  over  several  grid  points.  The  flux- 

corrected  transport  (FCT)  algorithm  (see  Reference  10) 

circumvents  this  difficulty  by  a  clever  computational 

scheme.  The  application  of  this  method  to  the  present 

problem  would  provide  benefits  in  dealing  with  shock  waves 

and  flow  perturbations  induced  by  area  change. 

The  grid  structure  of  the  present  problem  could  be 

streamline  by  constructing  a  truly  adaptive  grid.  This 

would  involve  transforming  the  governing  equations  with  a 

time-dependent  grid  transformation  equation  which 
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adequately  describes  the  grid  evolution  to  be  expected  for 
the  present  flow  field.  An  algorithm  to  cluster  grid 
points  at  locations  of  large  change  could  be  included  here. 
These  grid  calculations  would  require  the  solution  of  an 
additional  partial  differential  equation,  assuring  the 
continuity  of  the  transformation.  These  concepts  are 
discussed  in  References  11  and  12. 

VI.  RECOMMENDATIONS 

The  internal  flow  fields  of  the  light  gas  gun  are  of 
interest  from  both  operational  and  physical  viewpoints. 
It  is  recommended  that  the  baseline  performance 
prediction  model  discussed  in  Reference  5  be  developed 
to  its  full  potential  by  comprehensive  comparisons  with 
experimental  data  for  many  test  conditions.  Further, 
detailed  studies  of  the  flow  fields  within  the  various 
launcher  chambers,  in  the  manner  suggested  in  this  report, 
should  be  continued. 
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FIGURE  1.  TYPICAL  ARRABGEMEHT  OF  TOO-STAGE  LIGHT  GAS  GUM 
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PRESSURE-VOLUME  RELATIONS 
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Pc  ■  Coabuation  Chamber  Pressure 
N^  ■  Propellent  Weight  Burned 
Xp  ■  Piston  Velocity 
Xp  »  Piston  Displacement 
X,  »  Sabot  Velocity 
X8  ■  Sabot  Displacement 
Pjj  ■  Pump  Tube  Pressure 
P^  *  Launch  Tube  Pressure 
Y*  ■  Gas  Specific  Heat  Ratio 


TABLE  1.  BASELINE  PERFORMANCE  PREDICTION  EQUATIONS 
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